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ABSTRACT 


A study is made to improve the predictions of Reynolds stresses in a back- 
ward facing step flows, through modifications of the pressure-strain correla- 
tion. The mean-strain term of the pressure-strain correlation is formulated 
only in terms of non-isotropic turbulence in order to take the severe non- 
isotropic effect caused by a separating flow. This model is compared with 
other existing models and results are verified with existing experimental 
results. 
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NOMENCLATURE 


G 

Ga 

H 

B* 

k 

P 

Re 

S{j 

u 

U 

v 

x,y 


pressure strain correlation constants 
generation rate of turbulent kinetic energy 
generation rate of Reynolds stresses < U{Uj > 
step height 

secondary generation rate of Reynolds stresses 

< U{Uj > 

turbulent kinetic energy 
fluctuating pressure 

Reynolds number based on step height, H 
Surface integral 

fluctuating velocity in x direction 
mean velocity in x direction 
fluctuating velocity in y direction 
Cartesian coordinates 


Greek Symbols 

a, (3 

Sij 

e 

0ij,l > 

TT m * 

U V 

P 

u 


constants in Ily* 

Kronecker delta 

dissipation rate of turbulent kinetic energy 

pressure-strain correlations 

for Reynolds stresses 

fourth-order tensor 

density of fluid 

kinematic viscosity 


Subscripts / Superscripts 

i,j,k,l,m,n tensor notations 


Symbols 

0 


time-averaged values 



1 INTRODUCTION 


Separation and recirculation of flow is a common flow phenomena in many 
industrial problems. It associates itself with higher turbulence levels which 
not only render the flow greater analytical complexity but also result in 
highly augmented heat and momentum transfer. Therefore, mathematical 
models that can predict complex turbulent flows are needed. It has been 
commonly accepted that any isotropic turbulence models such as the k — 
e model et al. cannot take the non-isotropic effect into account in the 
evaluation of turbulence levels, whereas this effect is significant in complex 
turbulent flows. This is why the presents authors have been investigating 
to establishing a turbulence model based on the Reynolds stresses. The 
pressure-strain correlation term of the Reynolds-stress equations can be 
divided into two parts: the turbulence term and the mean strain term. The 
turbulence term controls the energy exchange between the normal stress 
components in relatively smaller eddies and the latter does the same in 
larger eddies. Our main focus is the evaluation of stresses in large eddies 
through the separating and reattaching flow which suffers strong curvature 
effect from the recirculation of flow. 

In improving the computations of the Reynolds stresses, the authors 
discovered that the closure models of the mean strain term of the pressure- 
strain correlation significantly affects the results in the separating and reat- 
taching flows. 

The model for the turbulence term in the pressure-strain correlation was 
developed by Rotta [1] and it has given effectively reasonable results with 
the anisotropy form. However, the existing models for mean strain term 
of the pressure-strain correlation proposed by several researchers [2-4] do 
not reach to a compromizing level in the computation of such a separating 
flow with recirculation. Thus, the authors have established a tensor that 
does not contain any terms which have isotropic turbulence but have only 
quantities with non-isotropic turbulence. The new model was compared 
with the existing models mentioned above. All models were then compared 
with experimental results of Chandrsuda and Bradshaw [5] with Reynolds 
number of 10 5 . 

Since constants of the existing models were determined by comparing 
with plane homogeneous shear flows [4], modifications were made for re- 
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circulating and separating flows through analytic comparison with experi- 
mental values , and computer optimazations. 


2 MATHEMATICAL MODELS 

2.1 Governing Equations 

The set of differential equations governing the transport of the kinematic 
Reynolds stresses ,< u t Uj > can be written in the form : 


d_ 

dxj 


( Uj < UiUj >) = — 


< UjUk > 


dUj 

dx k 


dU, 

< UiUk > q — 

dx k 


dui dit j 

- 2 v < - -- 1 > 


-v < 


du^dui > + < du^dui_ \ + < p ( dv^ 


dx{ dx 


dxj dxk 


dx k dx k 
dv,i \ 

+ 


p \dxi dx 


> 


_d_ 


< UiP > < UjP > 

< jU k ^ H ®jk “I - ®ik 


( d duk duk \ . 

> + + < ”'fc- >J (1) 

To close eq. (1) in terms of mean velocities and the Reynolds stresses, 
the turbulence quantities on the right-hand side of eq. (1) must be repre- 
sented by empirical functions of mean velocities, the Reynolds stresses and 
their derivatives. 

In this study, attention was focused on the approximation of the fourth 
group of terms of the right-hand side of eq. (1), namely the “pressure-strain 
correlation”. The rate of the “pressure-strain correlation” is the redistribu- 
tive effect ; thus, this correlation mainly controls the energy distribution. 

The authors have tested several models of the pressure-strain correla tion 
and discovered that it was necessary to remodel this correlation for the 
computations of flows which are accompanied by flow recirculation. Some 
existing models are summarized in the next subsection. 


2.2 Summary of Existing Models 

Generally, the pressure-strain correlation is divided into two parts: one 
with turbulence effect only and the other with both turbulence and the 
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mean strain rates. The first part with the turbulence effect was formulated 
by Rotta [1] in the form of anisotropy which is represented by the symbol 
and the last term has been proposed in several ways as follows : 


Model 1 : Naot et al. [2] 


p ( dui dui 

< ~ I J 


p \dxj ' dxj > ~ C ^ (GtJ 3 S,jG) + ^ 


Model 2 : Launder et al. [4] 


( 2 ) 


< - 


p ( dui du 


+ 


p \dxj ’ dxi 


>= - 


(C<t> 2 + 3) 


(3Q^-2) /c (dUi + du; 


11 

(8Q 2 - 2) 


(Ga - iS i3 G ) 


55 


dxi 


dxi 


[Hi j ^ijG) T i (3) 


where 



G = — < U{Uj 

dUi 

> dxj 

(4) 

Gij = — 

( du j 

< UiUk > Q - 

\ dx k 

dUA 

h<W> SxJ 

(5) 

H» = ~ ( 

( ou k 

< UiUk > a " 
V ° x j 

OUl, ' 

[- < UjUk > ^ 

| (6) 


C <f>\ ^ u i u j > 

(7) 


The first model was developed by removing the isotropic constraint from 
the double-velocity, two-point correlation tensor. The second model was 
obtained by approximating the pressure-strain correlation by an arbitrary 
fourth-order tensor, where the single-point correlation was used. 

The value of C ( p 1 was set equal to 1.5, and through computer opti- 
mizations, Cfc was recommended to be 0.7 and 0.4 for model of Naot and 
Launder respectively. 
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2.3 Formulation of A New Model 


Following Chou [6], a Poisson equation for the fluctuating pressure, p, can 
be obtained by taking the divergence of the equation for the fluctuating ve- 
locity U{, which can be re-expressed to write the pressure-strain correlation 
as 


p diii 

P dxj 


>— Sij+ 


47T 



> +2 



< 



> 


dvol 


x - y 


(8) 


where the first and second terms in the integral are denoted as 4>i^\ and 
2 respectively, and S v is a surface integral which is negligible away from 
a solid wall. 

Following Rotta’s proposal [1], 4>i h \ was formulated same as one given 
by eq. (7) with C* = 1.5. Furthermore, ^2 may be approximated as 

dUi 

^ (0) 

where n™‘ is a fourth-order tensor which satisfies the kinematic constraints 
of 


njy' = nj” 1 = nj (io) 

n^’ = o (11) 

n™‘ = 2 < u m ui > (12) 

Equation (12) suggests that II”" can be approximated by a linear combi- 
nation of the Reynolds stresses. On neglecting isotropy effects, the tensor 
can be expressed as : 

Ejj ' “ QiS t lj UiUj ^ UiU[ 
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f ^ tlm^j ^ ~k^ij ^ 'U'Tn'U'l ^ ) ”1” ^<f>2 til'll j ^ 


(13) 


where the condition of eq. (10) is already imposed. The application of cqs. 
(11), and (12) enables two of these constants to be expressed in terms of 
the third, namely Cfo and can be written as 

o — C* 2 + 1 




The final form is given as 


P ( duA 

p \dxj dxi ) 


>= -(2C^2 + 1 )(Gij - ^ SijG ) 


-2C t ,(H tj - \s xi G) - (14C * 3 1+2 b , ) G (14) 

The coefficient can be determined by simplifying the Reynolds 
stresses described by eq. (1) into an algebraic equation which was orig- 
inally developed by Rodi [7]. 

If the convection- diffusion string of the < upij > is set to be equal to 
that of the turbulence kinetic energy, k, it follows that, 


d 

dxj 


( Uj < UiUj >) — Dij = 


< U{Uj > 

k 



D k 


(15) 


where D,y and D k denote, respectively, the diffusion rates of < UiU 3 > and 
k. After some manipulation, eq. (15) can be re-arranged into the following 
form, 


UiUj > 



1 

C<f> J 6 



2C , 0 2 ((?,j 4- \ SjjG) + 2C^, 2 ( + jSjjG) 

C^e 


(16) 


5 


By assuming that the dissipative action is isotropic, 
parabolic, eq. (16) can be further reduced to : 


and the flow is 


{ <UU> ~ 2 3 



(6cy a + 1) 
C* 


(17) 


( 2 

< vv > — - 

V 3 



(6C*0 2 + |) 

C* 


( 18 ) 


(< ww > — 




(2 Cfo + I) 

C<t> i 


(19) 


The above set of relations are compared with typical experimental data 
for several parabolic flows, and the coefficient C'^ 2 is found as; 


C 


<t> 2 


1 
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Further investigations through comparison with recirculating flows, it was 
concluded that 


is the optimum value for prediction of separating and recirculating flows. 

3 NUMERICAL METHOD 

3.1 Numerical Procedure 

Formulation and discretization of all transport equations were performed by 
using the conventional control- volume approach of Patankar [8], by breaking 
each equation into diffusion, convection, and source terms. The system of 
equations were made linear so that they could be solved iteratively by the 
tridiagonal matrix algorithm. The convection-diffusion string was evaluated 
through the fourth-order approximation of the exponential form [9]. After 
several numerical tests, it was observed that the variation in results of 
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skin friction coefficient, pressure coefficient and mean velocity lie within 
2% when the grid is changed from 52 x 52 to 62 x 62. Therefore the grid 
independent state was assumed to be attained with the grid of 62 x 62 [10]. 

Computations are performed by using a variable grid mesh with the 
grid expanding linearly at the rate of 2% in the axial direction and at 3% 
in the transverse direction from the step wall. The iterative procedure is 
terminated when the maximum value of the relative residual sources of U, 
V, and mass balance falls below 1%, which is often achieved with 1000 
iterations. 

The complete process of solving the momentum, turbulent kinetic en- 
ergy and the Reynolds stresses takes about 60 minutes of CPU time on a. 
UNIVAC 1100 computer. 

3.2 Boundary Conditions 

At the inflow, outflow, and symmetry axis, conventional conditions are 
employed; prescribed values are provided at the inflow region, and zero 
gradients of the variable are used at the outflow region and the symmetry 
axis. 

At the wall boundary the wall function treatment is employed for the 
momentum and the turbulent kinetic energy. The dissipation rate is evalu- 
ated under the local equilibrium condition. Details of these conditions are 
given in [11]. The wall boundary conditions of the Reynolds stresses are as 
follow: 


- < uv >= -0.238k + (20) 

p ax 

<uu>=l.2\Ak (21) 

< vv >= 0.238 k (22) 

The above equations are used for the grids adjacent to a solid wall. 
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4 DISCUSSION OF RESULTS 

Figure 1 depicts the flow geometry and numerical grid used in the present 
computations. The height of the channel is magnified five times for better 
visualization. 

Figure 2 shows the velocity variations inside the solution domain. It 
can be seen that shear flow separates at the edge of the step and forms a. 
separating shear layer which causes a part of the fluid to recirculate in the 
region confined by a step and the bottom wall, and the rest to continue to 
flow downstream contributing to the stabilization of the flow through the 
formation of a boundary layer in the redeveloping region. 

Figures 3 and 4 compare the mean velocity profiles predicted by all three 
models at two different streamwise locations. It can be seen that the mean 
velocity is fairly independent of the choice of pressure-strain correlation. 

Figures 5 and 6 show the predicted shear stress. Generally all three 
models overpredicted the levels of the shear stress but the new model im- 
proved the prediction in the separating shear layer giving levels closer to 
the experimental data , whereas models by Naot et al. [2] and by Launder 
et al. [4] give levels of the Reynolds shear stress nearly zero. 

Figures 7 and 8 compare the < uu > profile. Although both models by 
Naot et al. and by Launder et al. give reasonable results, the agreement 
with the measured values in the separating shear layer is not quite satisfac- 
tory. The new model predicts well in this region but underpredicts at the 
wall. 

Figures 9 and 10 show the < vv > profile. The new model is far more 
superior than the other two models in predicting this Reynolds stress com- 
ponent in the separating shear layer, but does suffer from overprediction 
near the wall. 


5 CONCLUSIONS 

From the above results, it is clear that the development of a new pressure- 
strain correlation is essential to improve the prediction of the Reynolds 
stresses particularly in the separating shear layer. A near wall correction 
might be helpful to improve the prediction near the solid wall. 
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As mentioned before, the constants C < p s and C ^ might need further 
modifications to improve the predictions of Reynolds stresses. 
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FIGURE CAPTIONS 


Figure 1 Typical grid system used (5:1 magnification radially). 
Figure 2 Velocity vectors inside the solution domain. 

Figure 3 Mean velocity profile at fj = 5.4 

Figure 4 Mean velocity profile at fj = 6.4 

Figure 5 < uv > profile at ^ = 8.4 

Figure 6 < uv > profile at fj = 10.3 

Figure 7 < uu > profile at fj = 8.4 

Figure 8 < uu > profile at jj = 10.3 

Figure 9 < vv > profile at = 8.4 

Figure 10 < vv > profile at ^ = 10.3 
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